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Abstract — In  a  large  system,  such  as  a  water,  gas,  or 
electrical  distribution  system,  degraded  performance  due 
to  failures  of  components  can  be  modeled  as  a  set  of 
discrete  states  interconnected  by  edges  with  weights  that 
represent  conditional  probabilities.  To  establish  such  a  model, 
we  compute  the  conditional  probabilities  with  multi-valued 
decision  diagrams  (MDDs).  Since  a  typical  decision  diagram 
is  large,  the  computation  time  is  also  large.  In  this  paper,  we 
propose  an  edge- valued  MDD  (EVMDD)  based  method  to 
avoid  unnecessary  path  traversals.  The  proposed  method  is 
a  hybrid  method  of  a  path  traversal  method  and  a  bottom- 
up  method  that  visits  each  node  only  once.  By  effectively 
combining  both  methods,  we  achieve  a  speed-up  of  the 
analysis  by  about  2.3  times  for  large  systems  compared  to 
an  existing  method. 

Key  words  -  lYl  ulti -state  systems  with  multi-state  compo¬ 
nents;  structure  functions;  system  analysis  based  on  decision 
diagrams;  system  analysis  using  conditional  probabilities; 
EVMDDs. 

I.  Introduction 

Modem  various  systems,  such  as  computer  server  sys¬ 
tems,  telecommunication  systems,  water,  gas,  and  electri¬ 
cal  power  distribution  systems,  usually  achieve  tolerance 
to  faults  by  redundancy.  Even  if  a  fault  occurs,  these 
systems  still  keep  on  working  with  an  acceptable  or 
degraded  performance  level.  Thus,  modeling  such  systems 
by  two  states  (working  and  failure)  is  inadequate.  They 
are  often  modeled  by  multi-state  systems  with  multi¬ 
state  components  (or  simply  multi-state  systems)  in  which 
levels  of  performance,  reliability,  safety,  efficiency,  power 
consumption,  etc.  are  represented  by  more  than  two 
states  [2],  [3],  [12],  [14]-[  16]. 

To  design  fault  tolerant  systems,  intensive  analysis  of 
multi-state  systems  using  various  assessment  measures  for 
identifying  critical  components  and  system  weaknesses 
is  indispensable  [12].  Among  them,  assessing  the  prob¬ 
ability  of  each  state  of  a  multi-state  system  is  essential 
to  the  design  of  a  reliable  fault  tolerant  system  [14], 
[15].  Since  this  is  very  time-consuming,  many  methods 
to  shorten  analysis  time  have  been  proposed  [1],  [2], 
[4],  [6],  [10],  [13]— [15].  Especially,  methods  based  on 
decision  diagrams  have  attracted  much  attention  as  fast 
analysis  methods.  Unlike  conventional  graphical  models, 
such  as  Bayesian  networks,  that  represent  cause-and-effect 
relationships  among  components  and  a  system,  decision 
diagrams  represent  structure  of  a  system  as  a  function. 


Thus,  information  to  be  manipulated  is  reduced,  and  it 
results  in  faster  analysis. 

In  these  analysis  methods  of  multi-state  systems,  it  is 
assumed  that  states  of  components  in  the  systems  occur 
independently  of  other  components.  This  reduces  time 
complexity  [2],  [3]  of  the  analysis.  However,  there  exist 
systems,  such  as  medical  systems  [7],  in  which  states  of 
components  occur  depending  on  other  components.  Such 
systems  cannot  be  accurately  analyzed  by  the  methods 
assuming  only  independent  events  of  components.  We  may 
overlook  fatal  flaws  in  systems  by  using  such  analysis 
methods.  This  should  be  avoided  especially  in  safety- 
critical  systems  such  as  flight  control  systems  and  nuclear 
power  plant  monitoring  systems.  This  paper  shows  fast  and 
accurate  analysis  methods  for  multi-state  systems  having 
dependencies  among  components. 

In  [7],  a  multiple-valued  decision  diagram  (MDD)  based 
analysis  method  for  such  systems  has  been  proposed.  It 
can  analyze  systems  quickly  and  accurately  by  traversing 
paths  in  an  MDD  even  if  systems  have  dependencies 
among  components.  However,  systems  do  not  always 
have  dependencies  among  all  components.  There  can 
be  systems  in  which  only  a  part  of  components  has 
dependencies.  Thus,  in  this  paper,  we  propose  a  faster 
analysis  method  for  those  systems.  The  proposed  method 
is  based  on  an  edge-valued  MDD  (EVMDD)  [9],  [10]. 
It  avoids  unnecessary  path  traversals  of  EVMDDs  by 
applying  a  method  intended  for  only  independent  events  to 
components  without  dependencies,  and  analyzes  systems 
more  quickly  and  accurately. 

This  paper  is  organized  as  follows:  Section  II  defines 
multi-state  systems,  EVMDDs,  and  some  basic  terminolo¬ 
gies  related  to  probability.  Section  III  shows  the  analysis 
method  based  on  path  traversals  in  an  EVMDD,  and  in 
Section  IV,  we  propose  a  hybrid  method  to  avoid  unnec¬ 
essary  path  traversals.  Experimental  results  are  shown  in 
Section  V. 


II.  Preliminaries 

This  section  defines  multi-state  systems,  structure  func¬ 
tions,  EVMDDs  to  represent  structure  functions,  and  basic 
terminologies  related  to  probability. 
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Figure  1 .  Multi-state  system  for  an  electrical  power  distribution  system 
and  its  structure  function  [11], 


Figure  2.  MDD  for  the  structure  function  [11], 


A.  Multi-State  Systems  and  Structure  Functions 

Definition  1:  A  multi-state  system  is  a  model  of  a 
system  that  represents,  as  states,  a  capability,  such  as 
performance,  capacity,  or  reliability.  There  are  usually 
more  than  two  states,  and  so  a  multiple-valued  analysis 
is  required.  When  components  in  a  system  are  modeled 
as  well,  it  is  called  a  multi-state  system  with  multi-state 
components.  In  this  paper,  it  is  simply  called  a  multi-state 
system. 

Definition  2:  A  state  of  a  multi-state  system  depends 
only  on  states  of  components  in  the  system.  A  system  with 
n  components  can  be  considered  as  a  multi-valued  function 
f[x  i,X2,  . . . ,x„)  :  R\  x  i?2  x  . . .  x  R„  — >  M,  where  each  x; 
represents  a  component  with  r,  states,  If  =  {().  I ....  .r, 

1}  is  a  set  of  the  states,  and  M  =  {0, 1, . . .  ,m  —  1}  is  a 
set  of  the  m  system  states.  This  multi-valued  function  is 
called  a  structure  function  of  the  multi-state  system. 

Example  1:  Fig.  1(a)  shows  a  multi-state  system  for 
an  electrical  power  distribution  system.  In  this  figure,  the 
power  plants  xi,X2,x$  and  the  transformer  X4  have  three 
states  which  correspond  to  supply  levels:  0  (breakdown), 
1  (partially  supply),  and  2  (fully  supply).  And,  the  system 
has  six  states  which  correspond  to  the  percentage  of  area 
of  a  town  that  is  blacked  out:  0  (complete  blackout),  1 
(90%  blackout),  2  (60%  blackout),  3  (30%  blackout),  4 
(10%  blackout),  and  5  (0%  blackout). 

In  this  way,  by  assigning  a  value  to  each  state,  we 
obtain  the  6-valued  structure  function  /  shown  in  Fig.  1(b). 
Note  that  Fig.  1(b)  shows  a  part  of  the  34  =  81  entry 
table  since  it  is  too  large  to  be  included  in  its  entirety. 

(End  of  Example) 

B.  Edge-Valued  Multi-Valued  Decision  Diagrams 

Definition  3:  A  multi-valued  decision  diagram 
(MDD)  is  a  rooted  directed  acyclic  graph  representing 
a  multi-valued  function  /.  The  MDD  is  obtained  by 
repeatedly  applying  the  Shannon  expansion  to  the 
multi-valued  function  [5],  It  consists  of  non-terminal 
nodes  representing  sub-functions  obtained  from  /  by 
assigning  values  to  certain  variables.  It  also  has  terminal 
nodes  representing  sub-functions  obtained  from  /  by 
assigning  values  to  all  variables.  These  nodes  are  labeled 
by  function  values.  Each  non-terminal  node  has  multiple 


Figure  3.  EVMDD  for  the  structure  function  [11]. 


outgoing  edges  that  correspond  to  the  values  of  a  multi¬ 
valued  variable.  The  MDD  is  ordered;  i.e.,  the  order  of 
variables  along  any  path  from  the  root  node  to  a  terminal 
node  is  the  same.  The  MDD  is  reduced  by  applying  the 
following  two  reduction  rules: 

1)  Share  equivalent  sub-graphs. 

2)  Delete  non-terminal  nodes  whose  outgoing  edges  all 
point  to  the  same  node. 

When  an  MDD  represents  a  function  for  which  multi¬ 
valued  variables  have  different  domains,  it  is  a  heteroge¬ 
neous  MDD  [8],  In  the  following,  the  term  ‘MDD’  refers 
to  a  heterogeneous  MDD. 

Definition  4:  An  edge-valued  MDD  (EYMDD)  [9]  is 
an  extension  of  the  MDD,  and  represents  a  multi-valued 
function.  It  consists  of  one  terminal  node  representing  0 
and  non-terminal  nodes  with  edges  having  integer  weights; 
0-edges  always  have  zero  weights.  In  the  EVMDD,  the 
following  two  reduction  rules  are  applied: 

1)  Share  equivalent  sub-graphs. 

2)  Delete  non-terminal  nodes  whose  outgoing  edges  all 
point  to  the  same  node,  and  those  edge  weights  are 
0. 

In  the  EVMDD,  the  function  value  is  represented  as  a  sum 
of  weights  for  edges  traversed  from  the  root  node  to  the 
terminal  node. 

Example  2:  Fig.  2  and  Fig.  3  show  an  ordinary  MDD 
and  an  EVMDD  for  the  structure  function  of  Example  1. 


For  readability,  some  terminal  nodes  in  the  MDD  are  not 
combined.  (End  of  Example) 

C.  Conditional  Probabilities 
Definition  5:  A  conditional  probability,  denoted  by 
P(xi  =  Ci\xj  =  cf),  is  the  probability  that  an  event  jc,-  =  co¬ 
occurs  given  that  an  event  xj  =  cj  has  occurred.  The 
following  holds: 

n- 1 

X  p(xi  =  c< lx j  =  cj)=[-  (1) 

c,-=0 

Definition  6:  A  joint  probability,  denoted  by  P(x,  = 
Ci,Xj  =  Cj),  is  the  probability  that  both  of  an  event  Xj  =  Cj 
and  an  event  xj  =  cj  occur. 

The  conditional  probability  and  the  joint  probability 
can  be  extended  to  more  than  two  events.  The  following 
relations  hold  between  the  conditional  probability  and  the 
joint  probability: 

P(Xi  =  Ci,Xj  =  Cj)  =  P(Xi  =  Cj\Xj  —  Cj)  •  P(Xj  =  Cj).  (2) 
O'-' 

P(xi  =  a)  =  X  p(x>  =  C>’XJ  =  ci )’ 

Cj= 0 

o-1 

=  X  P(x‘  =  °i\Xi  =  Ci )  ■  ^(*7  =  C7'X3) 

c7=0 

Definition  7:  If  an  event  Xi  =  c,  does  not  directly  de¬ 
pend  on  an  event  Xj  =  c/  (i.e.,  the  events  are  conditionally 
independent),  then  the  following  relation  holds: 

P{xi  =  Ci,Xj  =  Cj\xk  =  ck ) 

=  P(Xi  =  Ci  \xk  =  Ck)-  P(Xj  =  Cj  \xk  =  ck) . 

Also, 

P(xt  =  Ci\Xj  =  C j,Xk  =  Ck)  =  P{Xi  =  Ci\xk  =  ck) . 

In  general,  thus,  the  joint  probability  is  computed  by: 

P(Xi  =Cl,X2=C2,  ...  1xn=c„) 

=  P{xi=cx\ii{xi))-P{x2  =  c2\n{x2))- 
■  ■■■P(xn=Cn  |II(jCn)),  (4) 

where  n(x,-)  is  a  set  of  variables  on  which  Xj  depends  and 
values  assigned  to  them. 

III.  Path  Traversal  Based  Analysis  Methods 

Definition  8:  The  probability  that  a  structure  function 
/  has  the  value  ^  is  denoted  by  Ps(f  =  s),  where  s  £ 
{0,  1,  ...  ,777  —  1}. 

An  analysis  of  multi-state  systems  is  to  solve  the 
following  problem: 

Problem  1:  Given  a  structure  function  /  of  a  multi-state 
system  and  the  conditional  probability  of  each  state  of 
each  component,  compute  the  probability  of  each  state  of 
the  multi-state  system  Ps{f  =  s). 

From  the  definition  of  the  joint  probability,  the  proba¬ 
bility  of  each  system  state  Ps(f  =  s)  can  be  computed  as 
follows: 

ps{f=s)=  X  P(f  =  s,xi,x2,...  ,X„),  (5) 

ceR " 


Figure  4.  Path  traversals  in  the  MDD. 


where  R"  is  a  set  of  value  assignments  to  all  the  variables 
Xi,  and  P(f  =  s,xi,X2,  •  ■  • ,  ,x„ )  is  a  joint  probability  of  /  = 
s  and  (xi,X2,-..,xn)  =c  =  (ci,C2,  ...,cn).  From  (2),  we 
have 

P(f  =  S,X  1  =Cl,*2  =  C2,...,Xn  =  C„)  = 

P(f  =  s\xi  =  Cl,X2  =  C2, . . .  ,Xn  =  Cn)  ■ 

P(x\  =  C\,X2  =  c2,...,xn  =  c„), 

where 

P(f  =  s \xi  =  c i,x2  =c2,...  ,x„  =  c„)  = 

0  (f(ci,c2,...,c„)^s) 

1  (/(ci,c2,...,c„)  =s) 

because  /  is  a  function  of  xi,x2, . . .  ,xn.  Therefore,  (5)  is 
obtained  by  a  sum  of  joint  probabilities  that  an  input  vector 
c  satisfies  f  =  s.  Since  each  joint  probability  can  be  com¬ 
puted  by  (4),  the  time  complexity  of  this  straightforward 
method  is  0(rn)  where  all  Xj  are  /--valued  variables. 

A.  Method  Based  on  Path  Traversals  in  MDD 

To  solve  Problem  1  more  efficiently,  a  method  that 
analyzes  systems  by  traversing  paths  in  an  MDD  has  been 
proposed  [7],  Since  the  number  of  paths  in  an  MDD 
becomes  smaller  than  r"  due  to  the  reduction  rule  for 
deleting  non-terminal  nodes,  this  analysis  method  is  faster. 

Example  3:  In  this  example,  assume  that  X4,xj,  and  x2 
depend  on  X},x2,  and  at,  respectively.  Then,  consider  path; 
in  Fig.  4.  Since  this  path  does  not  visit  a  node  for  at,  from 
the  reduction  rule  for  deleting  the  node,  we  compute  the 
following  using  conditional  probabilities  assigned  to  each 
edge  and  (1): 

2 

X  P(x 4  =  c4|a3  =  0)- 

C4— 0 

P(x i  =  0)  •  P(x2  =  0|at  =  0)  •  P(x3  =  0|a2  =  0), 

=  1  •  P( xi  —  0)  •  P(x2  =  0|ai  =  0)  •  P(x3  —  0|a2  =  0), 

=  P(x t  =  0,a2  =  0,a3  =  0). 

Similarly,  since  path2  does  not  visit  a  node  for  x3,  we 


Algorithm  1:  (Path  traversal  method  of  EVMDDs) _ 

l:path_EVMDD  (node  v,  vector  c,  edge  value  e,  probability  p )  { 

2:  if  (v  is  the  terminal  node)  { 

3:  prob_sys_state[e]  +=  p\ 

4:  return; 

5:  } 

6:  rj  =  domain  size  of  variable  x;  for  v; 

7:  for  (j  =  0;  j  <  r,-;  j++)  { 

8:  Cj  =  j; 

9:  p'  =  p  *  cond_prob(x,,  c,  1 

10:  path_EVMDD(j-th  edge  of  v,  c,  e  +  j-th  edge  value  of  v,  p') 

11:  } 

12:  Cj  =  undefined; 

m _ 


compute  the  following  using  (3): 

2 

X  p(-x3  =  c3^2  =  2)  •  P{x4  =  2\x3  =  c3)- 

C3— 0 

P(x  l  =2  )-P(x2  =  2\x\  =2), 

2 

=  X  pC*l  =  2’*2  =  2,x3  =  c3,x4  =  2), 

C3— 0 

=  P(x  i  =  2,X2  =  2,x4  =  2). 


(/'  =  .S)  =  P(Xj  =  Cj  I  X  ,=  fl  il  ■  P(Xj  =  dj) 

+  P(xj=  cj  x  j=  b  )  ■  P(Xi  =  bt) 

(a)  Xj  depends  on  x,\ 


P(xj  =  cj)  ■  { P(xt  =  at)  +  P(x,  =  b,)} 
(b)  xj  is  independent. 


In  this  way,  we  can  analyze  multi-state  systems  having 
any  dependencies  among  components  by  computing  a 
probability  on  each  path,  and  accumulating  it  at  each 
terminal  node.  (End  of  Example) 

B.  Method  Based  on  Path  Traversals  in  EVMDD 

In  this  subsection,  we  extend  the  MDD-based  method 
shown  in  the  previous  subsection  to  the  EVMDD-based 
one.  The  main  difference  is  only  the  handling  of  edge 
values.  While  the  MDD-based  method  accumulates  proba¬ 
bilities  on  paths  at  each  terminal  node,  the  EVMDD-based 
method  accumulates  probabilities  on  paths  in  each  sum  of 
edge  values.  Algorithm  1  shows  a  pseudo-code  for  the 
analysis  method  based  on  path  traversals  in  an  EVMDD. 

This  algorithm  requires  the  root  node  of  an  EVMDD, 
the  empty  vector  in  which  all  elements  are  undefined, 
the  edge  value  for  the  root  node,  and  the  probability 
for  the  root  node  1.0  as  input  parameters,  and  then  it 
traverses  all  paths  from  the  root  node  to  the  terminal  node 
recursively.  In  the  3rd  line,  the  array  prob_sys_state  stores 
a  probability  for  each  system  state.  In  the  9th  line,  the 
procedure  cond_prob(x,-,  c,  n(x,))  computes  a  conditional 
probability  P(xj  =  Ci|n(x,))  along  with  conditional  proba¬ 
bilities  for  skipped  parent  nodes,  as  shown  in  Example  3. 
Note  that  a  variable  x,-  depending  on  other  variables  ITfx,) 
has  to  be  below  those  variables  n(x,)  in  the  variable  order 
of  EVMDD  to  compute  conditional  probabilities. 


Figure  5.  Computation  of  conditional  probabilities  in  EVMDDs. 
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Figure  6.  Partition  of  EVMDD. 


is  shared  due  to  the  reduction  rule  for  sharing  equivalent 
sub-graphs.  As  shown  in  Fig.  5(a),  when  Xj  depends  on 
x;,  a  conditional  probability  has  to  be  computed  in  each 
path  since 

P(xj  =  Cj \xt  =  dj)  ■  P(xj  =  di)  +  P{xj  =  Cj\xi  =  bj)  ■  P(xj  =  bi) 

7 L  P{xj  =  Cj\xj  =  dj)  •  {P(xi  =  dj)  +P(xi  =  bi)}. 

On  the  other  hand,  when  xj  is  independent  of  others  as 
shown  in  Fig.  5(b),  probabilities  obtained  from  parent 
nodes  can  be  tallied  at  each  node  since 


IV.  Hybrid  Analysis  Method  for  Systems 
Partially  Having  Dependent  Components 
Although  Algorithm  1  can  analyze  multi-state  systems 
having  any  dependencies  among  components,  systems  do 
not  always  have  dependencies  among  all  components. 
There  can  be  systems  in  which  only  some  of  the  com¬ 
ponents  have  dependencies.  For  analysis  of  such  systems, 
Algorithm  1,  which  is  based  on  path  traversals  in  an 
EVMDD  is  not  very  efficient.  This  is  because  the  path 
traversal  based  methods  can  visit  a  node  many  times  that 


P(xj  =  Cj)-P{xi  =  di)+P{xj  =  Cj)-P(xi  =  bi) 

=  P(x j  =  Cj)  ■  {P(xj  =  dj)+P(xi  =  bi)}. 

Thus,  we  can  apply  a  method  that  visits  each  node 
only  once  to  components  without  dependencies  to  avoid 
unnecessary  path  traversals  of  EVMDDs. 

Suppose  that  an  EVMDD  can  be  partitioned  into  two 
parts:  the  upper  and  the  lower  parts,  as  shown  in  Fig.  6, 
where  the  upper  part  has  components  that  depend  on 


Figure  7.  Analysis  using  the  hybrid  method. 


others,  and  the  lower  has  only  independent  components*. 
Then,  we  apply  the  path  traversal  based  method  to  the 
upper  part,  and  apply  the  node  visits  based  method  [10] 
to  the  lower  part  in  a  bottom-up  manner.  The  method  [10] 
visits  each  node  only  once  in  a  bottom-up  manner.  And  at 
each  node,  it  computes  probabilities  for  its  sub-function 
by  merging  probabilities  obtained  at  child  nodes. 

We  begin  by  applying  the  node  visit  based  method  to 
the  lower  part.  Although  in  [10],  the  method  is  applied  to 
all  components  from  the  terminal  node  to  the  root  node, 
in  the  hybrid  method,  it  is  applied  only  to  independent 
components.  Then,  we  apply  the  path  traversal  based 
method  to  the  upper  part.  When  a  path  traversal  reaches  a 
node  v  in  the  lower  part,  the  probability  p  on  the  path  is 
multiplied  by  probabilities  Ps  (  f\,  =  ,v)  for  the  sub-function 
of  v'  that  are  obtained  by  the  node  visit  based  method. 
The  products  p-Ps(fv  =  s)  are  accumulated  as  probabilities 
Ps(f  =  s  +  e)  of  system  states  s  +  e,  where  e  is  a  sum  of 
edge  values  the  path  traversed. 

Example  4:  In  this  example,  assume  that  x2  depends  on 
X] ,  and  there  are  no  other  dependencies.  Then,  as  shown  in 
Fig.  7,  nodes  for  x\  and  x2  are  in  the  upper  part,  and  nodes 
for  X3  and  X4  are  in  the  lower  part.  First,  probabilities  at 
each  node  in  the  lower  part  are  computed  in  a  bottom-up 
manner  [10].  Next,  probabilities  on  paths  from  the  root 
node  to  nodes  in  the  lower  part  are  computed.  Consider 
pathi  in  Fig.  7.  This  path  reaches  node  vi,  and  thus  the 
probability  pi  =  P(x 2  =  0|xi  =  0)  -P(x \  =  0)  on  the  path 
is  multiplied  by  probabilities  Ps(fVl  =  0)  and  Ps(fVl  =  1) 
at  vi,  respectively.  The  products  pi  ■  Ps(fVi  =  0)  and  pi  ■ 
Ps(fvi  =  1)  are  accumulated  as  probabilities  Ps(f  =  0)  and 
Ps(f  =  1),  respectively,  since  the  sum  of  edge  values  e  for 
the  path  is  0.  Similarly,  since  path2  reaches  node  V2,  the 
products  p2  ■  Ps (fv2  =  0),  P2-Ps{fv2  =  3),  and  p2  ■  Ps(fVl  = 
5)  are  accumulated  as  probabilities  Ps(f  =  0),  Ps(f  =  3), 
and  Ps(f  =  5),  respectively. 

In  this  way,  we  analyze  systems  having  partially  de¬ 
pendent  components  efficiently  by  combining  the  path 
traversal  and  node  visit  based  methods.  (End  of  Example) 

*If  the  lower  part  has  dependent  components,  the  variable  order  is 
changed  so  that  the  lower  part  has  only  independent  components. 


Computation  on  each  path  in  the  path  traversal  based 
method  is  simpler  than  computation  at  each  node  in 
the  node  visit  based  method.  Thus,  the  path  traversal 
based  method  is  faster  when  the  number  of  paths  is  not 
much  larger  than  that  of  nodes  (i.e.,  when  the  number  of 
shared  nodes  is  small).  As  the  shape  of  the  EVMDD  in 
Fig.  6  shows,  nodes  in  the  upper  part  are  hardly  shared 
while  nodes  in  the  lower  part  tend  to  be  shared.  Thus, 
by  applying  the  path  traversal  and  the  node  visit  based 
methods  to  the  upper  and  the  lower  parts,  respectively, 
we  can  analyze  systems  efficiently.  In  the  experiment  in 
Section  V,  we  partitioned  an  EVMDD  with  n  components 
into  the  upper  part  with  |~|n]  components  and  the  lower 
part  with  components  because  it  produces  the  best 
performance  for  a  large  system. 

V.  Experimental  Results 

To  show  the  efficiency  of  the  proposed  hybrid  method, 
we  compare  it  with  various  existing  methods  using  the 
same  structure  functions  as  [10].  This  is  because  unfor¬ 
tunately,  benchmark  structure  functions  of  large  enough 
multi-state  systems  are  unavailable.  In  this  experiment,  we 
analyze  the  multi-state  systems  given  by  those  structure 
functions,  assuming  that  only  a  component  x2  depends  on 
another  component  x2-  The  methods  are  implemented  on 
our  private  EVMDD  package,  and  run  on  the  following 
computer  environment:  CPU:  Intel  Core2  Quad  Q6600 
2.4GHz,  memory:  4GB,  OS:  CentOS  5.7,  and  C-compiler: 
gcc  -02  (version  4.1.2).  Table  I  shows  the  experimental 
results  for  randomly  generated  m-state  systems  with  n  3- 
state  components. 

From  this  table,  we  can  see  that  the  methods  based 
on  path  traversals  in  MDDs  and  EVMDDs  are  much 
faster  than  the  naive  method  that  is  used  in  Bayesian 
networks  [16].  This  is  because  the  number  of  paths  in 
MDDs  and  EVMDDs  is  much  smaller  than  3”.  The 
method  based  on  path  traversals  in  EVMDDs  is  slightly 
slower  than  that  in  MDDs  because  of  the  overhead  to 
compute  edge  values. 

For  large  systems,  the  hybrid  method  using  both  path 
traversal  and  node  visit  based  methods  is  faster  than  the 
path  traversal  based  methods.  Especially,  as  a  difference 
between  the  number  of  paths  and  the  number  of  nodes  in 
an  EVMDD  becomes  larger,  the  hybrid  method  becomes 
faster.  In  this  experiment,  when  the  number  of  nodes  is 
smaller  than  ^  of  the  number  of  paths,  the  hybrid  method 
is  faster  than  the  path  traversal  based  methods.  On  the 
other  hand,  in  ordinary  MDDs,  the  number  of  nodes  is 
not  much  smaller  than  the  number  of  paths,  since  the 
number  of  shared  nodes  is  small.  Thus,  in  this  experiment, 
we  could  not  achieve  a  speed-up  by  applying  the  hybrid 
method  to  MDDs. ' 

From  the  above  observations,  we  can  conclude  that  the 
efficiency  of  the  hybrid  method  is  due  in  large  part  to  the 
compactness  of  EVMDDs,  and  it  is  the  method  that  can 
take  advantage  of  the  fact  that  EVMDDs  are  compact. 

^Results  about  those  are  omitted. 


Table  I 

Computation  times  for  analysis  of  w-state  systems  with  n  3-state  components. 


n 

m 

Number  of  nodes 

Number  of  paths 

Computation  time  (ju sec.) 

MDD 

EVMDD 

MDD 

EVMDD 

Naive  (5) 

Path  traversals  methods 

Ratio 

MDD  [7] 

EVMDD 

5 

3 

WMSM 

HI 

131% 

5 

10 

■a 

mg 

140% 

10 

3 

17 

15 

59,049 

29 

29 

4,030.00 

1.98 

2.01 

2.74 

139% 

10 

10 

77 

57 

135 

135 

4,019.00 

8.83 

9.12 

12.14 

137% 

10 

100 

599 

265 

999 

999 

4,011.00 

63.28 

66.35 

69.60 

10 

6,403 

4.031.00 

419.51 

437.76 

360.14 

86% 

15 

3 

32 

14,348,907 

59 

59 

1,356.000.00 

3.77 

3.91 

5.95 

158% 

15 

10 

14,348,907 

221 

221 

1,358.000.00 

14.84 

21.14 

151% 

15 

100 

14,348,907 

1,997 

1,997 

1,359.000.00 

132.65 

148.64 

119% 

15 

1,000 

3,363 

14,348,907 

1,360.000.00 

930.28 

92% 

15 

10,000 

11,474 

14,348,907 

1,360.000.00 

5,408.00 

15 

495,224 

62,759 

14,348,907 

1,359.000.00 

24,238.00 

43% 

n:  Number  of  3-state  components.  m:  Number  of  states  for  systems.  Ratio:  Hybrid  /  MDD  [7]  x  100  (%) 

Naive  (5):  The  method  that  directly  computes  (5). 

Very  short  computation  times  are  average  times  obtained  by  running  the  same  computation  1,000,000  times,  and  dividing  its  total  time  by  1,000,000. 


VI.  Conclusion  and  Comments 

This  paper  proposes  an  efficient  analysis  method  of 
multi-state  systems  in  which  states  of  some  components 
occur  depending  on  states  of  other  components.  The 
proposed  method  is  a  hybrid  method  using  both  path 
traversals  and  node  visits  in  an  EVMDD.  By  using  the 
node  visit  based  method,  we  can  avoid  unnecessary  path 
traversals  and  achieve  fast  analysis.  Experimental  results 
show  that  we  can  analyze  large  systems  quickly  and 
accurately  using  the  proposed  hybrid  method.  For  a  large 
system,  we  achieve  a  speed-up  of  2.3  times,  compared  to 
the  existing  path  traversal  based  method.  For  future  work, 
we  will  study  a  heuristic  method  to  find  a  better  partition 
of  an  EVMDD  in  Fig.  6  to  improve  performance  of  the 
hybrid  method. 
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